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Abstract 

We consider spatially homogeneous, anisotropic cosmological models in 5D whose line element can be written 
as dS^ — A{u,v)dudv — Bij{u,v)dx^dx-' , = 1,2,3), where u and v are light-like coordinates. In the case 
where Bij is diagonal, we construct three families of analytic solutions to the 5D vacuum field equations Rab = 
{A,B = 0,1,2,3,4). Among them, there is a family of self-similar homothetic solutions that contains, as a 
particular case, the so-called light-like Kasner universes. In this work we provide a detailed study of the different 
types of AD scenarios that can be embedded in such universes. For the sake of generality of the discussion, and 
applicability of the results, in our analysis we consider the two versions of non-compactified 5D relativity in vogue, 
viz., braneworld theory and induced matter theory. We find a great variety of cosmological models in 4D which 
are anisotropic versions of the FRW ones. We obtain models on the brane with a non-vanishing cosmological 
term A(4), which inflate a la de Sitter without satisfying the classical false- vacuum equation of state. Using the 
symmetry of the solutions, we construct a class of non-static vacuum solutions on the brane. We also develop 
static pancake- like distributions where the matter is concentrated in a thin surface (near z — 0), similar to those 
proposed by Zel'dovich for the shape of the first collapsed objects in an expanding anisotropic universe. The 
solutions discussed here can be applied in a variety of physical situations. 
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1 Introduction 



In recent years there has been an increased interest in theories that envision our spacetime as embedded in a universe 
with more than four large dimensions. There are several reasons that justify this interest, among them that extensions 
of four-dimensional general relativity to five and more dimensions seem to provide the best route to unification of 
gravity with the interactions of particle physics In 5-D there are two versions of relativity where the extra 

dimension is not assumed to be compactificd. These are membrane theory [5] and space-time-matter (or induced 
matter) theory [6j . They lead to a great variety of models both in the cosmological context and in the description of 
local self-gravitating objects (see, e.g., [7], [5]). Most of these models have been obtained in coordinates where the 
metric in 51? can be written a^ 

dS^ = g^^{xP, '4,)dx^'dx'' + e<^'^{xP, (1) 

in such a way that our AD spacetime can be recovered by going onto a hypersurface : = ipo — constant, which 
is orthogonal to the 5D unit vector 



A 



n 



A _ ^4 „ „A 



UAn"^ = e, (2) 



along the extra dimension, and g^^, can be interpreted as the metric of the spacetime. 

In this framework, the effective equations for gravity in AD are obtained from dimensional reduction of the 
Einstein field equations in 5D. The reduction is based on Campbell's theorem [10], [11] and consists in isolating 
the AD part of the relevant bD geometric quantities and use them to construct the AD Einstein tensor '^'Gq/j. The 
crucial result is that, even in the case where the energy-momentum tensor (EMT) in 5-D is zero, to an observer 
confined to making physical measurements in our ordinary spacetime, and not aware of the extra dimension, the 
spacetime is not empty but contains (effective) matter whose EMT, ^^-'Tq,^, is determined by the Einstein equations 
in AD, namely 

(4)G„0 = 8^ (4)r<,^ = {K^y^K^p ~ KlK^p) + ^g^p {KxpK^P - {Klf) - eE^p, (3) 



where Knu is the extrinsic curvature 



2 "^"'^ 2$ dij 

E^y is the projection of the bulk Weyl tensor ^^"^C abcd orthogonal to n"^, i.e., "parallel" to spacetime, viz., 

Eap 'C^ApBn n ^ --—^ + KapKp- e—^, (5) 

and = d^/dx°'. Before going on, it is worthwhile to emphasize that the above dimensional reduction of the 
field equations in bD is a standard technique that leads to the same effective matter content in AD, i.e., ^^"^Tap, 
regardless of whether the line element ([T]) is interpreted within the context of brane theory with Z2 symmetry |12] 
or space-time-matter (STM) theory lj5 . In this sense these two approaches to bD relativity are mathematically 
equivalent. However, they are different as regards physical interpretation and motivation [14 . In brane theory there 
is a singular hypersurface that defines spacetime, and the properties of matter in that hypersurface are, in general, 
not identical to the ones of induced matter calculated in STM from the effective EMT defined by ([3]). 

In the cosmological realm, nearly all models assume spatial homogeneity and isotropy, which means that the line 
element in bD is taken to be an extended version of the conventional Friedmann-Roberson- Walker (FRW) metric in 
AD, namely 

dS'^ = n^{t,i;)dt^ ~a'^{t,iP)-f,jdx'dx^ +e^^{t,i;)di;^, = 1,2,3, (6) 



Notation: = (x'^ , , , x^) are the coordinates in 4D and 1/1 is the coordinate along the extra dimension. We use spacetime 
signature — , — , — ), while e = ±1 allows for spacelike or timelike extra dimension, both of which are physically admissible for a detailed 
discussion see, e.g., [9]. 
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where 7^ is a maximally symmetric 3-dimensional metric, with curvature index fc = 0, ±1. In these coordinates the 
full integration of the vacuum Einstein field equations in 51? requires the specification of two additional assumptions. 
One of them is usually an assumption of geometric nature, e.g., that $ = 1, or n = 1. The second one is usually an 
equation of state for the matter quantities in 4D [TS], [TO] . 

Observations indicate that on large scales 100 Mpc) the universe is homogeneous and isotropic and well 
described by spatially-flat FRW cosmologies. However, there is no reason to expect that such features should hold at 
the early stages of the evolution of the universe. Rather, it is generally accepted that anisotropy could have played 
a significant role in the early universe and that it has been fading away in the course of cosmic evolution. In the 
framework of 4-dimensional spacetime, a prototype for anisotropic vacuum cosmologies is provided by the Kasner 
metric [17] , which mimics the behavior of more general solutions near the singularity during some finite periods of 
tim^. Various higher dimensional extensions of the vacuum Kasner model have been discussed in the literature [18j , 

m- 

In this paper we consider spatially homogeneous but anisotropic cosmological models whose metric in 5D has the 
form 

dS^ — A{u, v)dudv — Bij{u, v)dx^dx\ (7) 

where A and B are some functions of the "light-like" coordinates u and v. These metrics are different from ([6]) in 
various aspects: (i) they do not contain the time or extra dimension in an explicit way; (ii) the hypersurfaces of 
constant u or v are three-dimensional instead of 4Z3; (ii) a priori it is not clear how to define the 5D unit vector 
fi^ along the extra dimension, which in turn is needed for defining the spacetime sections and for constructing the 
appropriate projected quantities in 4D. 

Here, for the case where Bij is diagonal we construct three families of analytic solutions to the bD vacuum field 
equations Rab — {A,B — 0, 1,2,3,4). The simplest one is a family of self-similar solutionf0, which contains as 
a particular case the so-called light-like Kasner universes. From a physical point of view, self-similar homothetic 
models are interesting because they may serve as asymptotic regimes, i.e., near the initial cosmological singularity 
and at late times, for many homogeneous and inhomogcneous cosmological models |24] . The other two families of 
solutions are obtained under the assumption that some of the metric coefficients are separable functions of their 
arguments. 

In view of their potential relevance to the "similarity hypothesis" [24] , in this work we focus our attention to the 
family of self-similar 5Z> spacetimes mentioned above. The main question under study is what kind of AD scenarios 
can be embedded in such spacetimes. For the sake of generality of the discussion, and applicability of our results, 
in our analysis we consider both versions of non-compactified bD relativity, viz., induced-matter and brane theory. 
Unfortunately, the expressions obtained in AD as projections of the 5£) solutions are quite complex and cumbersome. 
Therefore, to obtain manageable mathematical expressions in 4D, in sections 3, 4 and 5 we simplify the algebra (but 
not the physics) by restricting our discussion to the subset of light-like Kasner solutions. 

Our solutions generalize a number of isotropic cosmological models and give back previous ones in the literature 
(see e.g., [25]- [27] and references therein). Although we are not discussing particular applications here, they could 
be useful in the study of generalizations of Mixmaster or Belinskii-Khalatnikov-Lifshitz oscillations in theories with 
a single extra dimension [21], [22], [28]- [30]. They could also be apphed to studying conjectures about isotropic Big 
Bang singularities in braneworlds [31]- [33]. Certain cosmological models, such as the cyclic universe model, also 
require an understanding of the behavior of Kasner-like solutions and are based on brane- type models [34j . 

The paper is organized as follows. In section 2 we derive our self-similar solution (the other two families of 
solutions are presented in the Appendix) and introduce a timelike coordinate t and a spacelike coordinate along 
the extra dimension. This is equivalent to introducing two additional degrees of freedom, which are expressed in 
terms of two functions of t and i/;. We will see that, as in the familiar FRW picture ([6]), these two functions can be 
related to the specific choice for embedding E^, in 5D and to the physics in AD. In section 3, within the context of 
STM we show that the light- like Kasner solutions generate a great variety of cosmological models in 4D, including 

^Without entering into technical details: extrapolating backwards in time towards the singularity, one finds an infinite number of 
alternating quasi-periodic Kasner-like epochs with different expansion rates I20I - I22I . 

^In the traditional interpretation of Sedov, Taub and Zel'dovich 1231 . self-similarity means that all dimensionlcss quantities in the 
theory can be expressed as functions only of a single similarity variable, which is some combination of the independent coordinates. In 
this way the field equations become a system of ordinary, instead of partial, differential equations. 
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the de Sitter, Milne and power-law FRW models. In section 4, within the context of the braneworld paradigm we 
find that they embed AD cosmological models with a non- vanishing cosniological term A(4') , which in principle can 
be either constant or time-dependent. In the case of constant A(4) the 3D space exponentially inflates, regardless 
of the specific embedding. We also show that by virtue of the symmetry {x,y,z) 4-> ?/;, they generate a class of 
non-static vacuum solutions on the brane. In section 5, once again using symmetry properties, we demonstrate that 
the Kasner-like metric ([22|) can be used to generate static pancake-like distributions of matter in AD, similar to 
those proposed by Zel'dovich for the shape of the first collapsed objects in an expanding anisotropic universe [35j . 
In section 6 we present a summary of our results. 



2 Cosmological models in 5D. Light-like coordinates 
2.1 Solving the field equations. Part I 

In this work we obtain three families of solutions to the field equations Rab = 0. However, to facilitate the 
discussion, in this subsection we present only one of them. Specifically, we present the family of solutions that we 
will use throughout the paper, which (as we will see in sections 3-5) may be interpreted or used as 5D embeddings 
for a number of AD universes. The derivation of the other two families of solutions, whose AD interpretation is not 
discussed here, is deferred to the Appendix (Solving the field equations. Part II). 

To simpfify the shape of the field equations let us momentarily denote Bn — e^^"'^\ B22 — e''^"'"' and S33 = 
gcr(«,f)_ From R,j~x = 0, Ryy = Rzz = we obtain the equations 

AXuv + Xu i<yv + + fJ-v) + ifJ-u + CTu) = 0, 

A^uv + ^J'uiK + '2^J'v + o-y) + ^J'v {o-u + K) = o, 

A(7uv + cTu {pv + '2ay + Xy) + ay {Xu + fJ-u) = 0. (8) 

Here the subscripts u, v indicate partial derivatives with respect to those arguments. The above equations show 
cyclic permutation symmetry, i.e., starting from any of them by means of the transformation A — *■ /i — > tr — > A we 
obtain the other two. 

• Self-similar solutions: First we solve the field equations under the assumption that the metric ([7]) possesses 
self-similar symmetry. This assumption is motivated by a number of studies suggesting that self-similar models 
play a significant role at asymptotic regimes [24] . From a mathematical point of view, it means that by a suitable 
transformation of coordinates all the dimensionless quantities can be put in a form where they are functions only of 
a single variable (say 23J. In our particular case, this implies that A = A(C), ^ = /^(C), c = ""(C), where is some 
function of u and v, viz., 

C = (9) 
With this assumption the first equation in ([8]) reduces to 

2^ + (A,+^.,+ac) + 2-^ = 0. (10) 

The assumed symmetry requires (Cuv / CuCv) to be some function of C, say Z((^) — {Cuv / CuCv) ■ Integrating we get 

Ac = 2a / ^(O^^C ^ 2a (Jj^, (11) 

where a is an arbitrary constant of integration. Similar equations, with new constants, e.g., (3 and 7, are obtained 
for fi and a by means of a cyclic transformation. What this means is that 

^ ^ ^ ^ k. C121 
2a 2/3 27 / ' ^ ' 



which upon integration yields 



e 



^ = €^^02^, e'^^Csf^, (13) 
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where Ci, C2, C3 are constants of integration. A single differential equation for f{() — f{u, v) can be easily obtained 
by substituting (|13p into any namely 

ffuv + (a - 1) fufv = 0. (14) 

Notation: Here and henceforth we denote 

a = a + /? + 7, 6 = + + 7^ - a - /? - 7, c = af] + + p-f, (15) 

where a, /3, 7 are arbitrary parameters. 
A simple integration gives 

f = [hiu) + g{v)f\ (16) 

where h{u) and g{v) are arbitrary functions of their arguments. Clearly, in the present case Bij are power-law type 
solutions of the similarity variable ^ = [h(u) + g(v)]. 

To simplify the discussion, and eliminate spurious degrees of freedom, we now make the coordinate transformation 
h(u) = ciu, g{v) = C2V, where Ci and C2 are constants. In these new coordinates 

A{u, v)dudv A{u, v)dudv, (17) 

where A{u,v) = [A(u,v)/hugv] with u and v expressed in terms of u, v. Then relabeling the coordinates (dropping 
the overbars) the metric becomes 

dS^ = A{u, v)dudv - Ci [cm + C2vf°'''' dx^ - C2 [cm + C2vf^''' dy^ - C3 [cm + C2vf'^l'' dz^. (18) 
For this metric the field equations — and R^y — yield the equations 

[cm + C2V] An + 2c CiA = 0, 

[cm + C2v] Ay + 2c C2A = 0, (19) 

which have a unique solution given by ^ = Co {cm + C2v) , where Co is a constant of integration. Now, it is 
easy to verify that = is identically satisfied. In summary, the final form of the self-similar solution is given bj{f| 

A={cm + c2v)-^'^''\, 622=^-''"/^ 633 =-4-^"/^ B,,=Q, i^j. (20) 

We note that this solution admits a homothetic Killing vector in 5D for any values of a, (3 and 7, namely, 

CcgAB = 2gAB, with = [t]qu, r]ov, (1 - a'qo/a)x, (1 - l3riQ/a)y, (1 - jriQ/a)z] , (21) 

where gAB is the metric (|20p . denotes the Lie derivative along the 5D vector and ryo = ci^ /{a^ ~ c). In addition, 
by setting one of the constants equal to zero, say C2 — 0, and making the coordinate transformation u~ du du, 
it reduces to 

dS^ = dudv -AvP^dx^- BvP^ dy^ - CvP^ dz^ , (22) 
where A, B, C are constants with the appropriate units, and pi,p2,P3 denote 

_ 2a{a + (3 + j) ^ 2(i{a + /3 + 7) ^27(0+^3 + 7) 

a2 + /32 + 72 ' ^2 ^2 + ^2 ^ ^2 ' P3 ^2 ^. ^2 ^2 ' K^-^) 

which satisfy the relation Yl^i=i (Pi ~ 1)^ = ^ '^'^2/ values oi a, f3 and 7. The metric ([^^ is usually called light-like 
Kasner solution. In this case the 5D homothetic vector is given by C*" = [u, v, (1 — pi/2)x, (1 — p2/2)y, (1 — p3/2)z]. 



The proportionality coefficients Co, Ci, C2, C3 can be set equal to unity without any loss of generality. 
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2.2 Introducing the timelike and "extra" coordinates 

In order to be able to apply the standard dimensional reduction ^ to metrics ([7]) one has to introduce coordinates 
that are adapted to the spacetime sections S^. With this aim we make the coordinate transformation 

u^F{t,ij), v^V{t,i^), (24) 

where t is assumed to be the timelike coordinate; the "extra" coordinate; F and V are, in principle, arbitrary 
differentiable functions of their arguments, except from the condition that the Jacobian of the transformation must 
be nonzero. 

With this transformation we obtain 

dudv = (FVdt^ + F'V'd^^) + (fV' + F'v) dtd^p, (25) 



where dots and primes denote derivatives with respect to t and ip, respectively. We can choose the coordinates t, ip 
in such a way that the 5D metric be diagonal. This requires 



F'V 

V = (26) 

F ^ 



As a consequence, the line element ([7]) becomes 



= A[t, ij)FVdt^ - B,j{t, ^p)dx'dx^ - Ait, ip)^^dip^, (27) 

F 

where A{t,ip) = A{F,V) and Bij{t,ifj) = Bij{F,V). A couple of points should be noticed here. Firstly, that the 
physical requirement goo > demands ip to be spacelike. Secondly, that the line element contains two arbitrary 
functions, which are not present in the original solution (|20p . The question is, why? Is this a mathematical, or gauge, 
artifact? 

The answer to this question is that the arbitrary functions in (I27|) arc not gauge artifacts. They reflect the 
physical reality that there are many ways of embedding a 41? spacetime in 51? while satisfying the field equations. 
If we choose some particular embedding we obtain a differential constraint connecting V and F, which allows us to 
obtain one of them in terms of the other, e.g., V in terms of F. Then, the remaining unknown function, e.g. F, 
can be determined from the physics in AD. As an illustration of the former assertion, let us consider two common 
embeddings that arise from the choice of the coordinate/reference system. 

Gaussian normal coordinate system: A popular choice in the literature is to use the five degrees of coordinate 
freedom to set g^^ — and (744 = —1. This is the so-called 'Gaussian normal coordinate system' based on E^. 
Consequently, in such coordinates V — {F/AF'^) and (pS)) becomes V ~ —{1/AF'). Now the condition (dV/d^p) = 
[dV'Idt] yields 

(ciF + C2V^)F" - 2cciF'2 0. (28) 
If C2 7^ 0, this equation gives V{t, %p) for any smooth function F{t, ip), and the metric ((27)) becomes 

dS^ = (^y^ dt^ - B.jit, tp)dx'dx^ - dip'^. (29) 
If C2 = 0, then is an equation for F. Integrating it we find 

F ^[l[t)+i;h(t)f''^''+^\ (30) 
where lit) and h{t) are arbitrary differentiable functions. 
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Synchronous reference system: The choice 500 = 1 is usual in cosmology: it corresponds to the so-called 
synchronous reference system where the coordinate t is the proper time at each point. Thus, setting V — (\IAF\ 
the line element (P7)) becomes 

dS^ = dt^ - B,j{t, iP)dx'dx^ - (^y^ diP"^. (31) 

In these coordinates (gH) reduces to V = -{F'/AF'^) and {dV/diJj) = (dV'/dt) yields 

[ciF + C2V) F - 2cciF2 ^ q. (32) 

Thus, for C2 = we get 

F = [M{i;) + tiV(V')]"'/("+'') , (33) 

where M and N are arbitrary differentiable functions of ip. For any other C2 7^ 0, we obtain V from ((32)l after 
choosing some smooth function F{t^ tjj). 

Thus, in principle the function V can be determined if we know F. At this point the question arises of whether 
we can single out the function F from "physical" considerations in AD. Further analysis of the field equations shows 
that if we assume an equation of state for the matter in AD, then we obtain an extra differential equation connecting 
V and F, which in addition to (|28|) or ([32|) . allows to express the solution (|27|) in terms of t and ■(/;. This is what is 
required for the 4+1 dimensional reduction of the 5D solutions. The general calculations are straightforward, but 
the equations are notational cumbersome in both STM and braneworld theory. On the other hand, ([50)1 and ([55)1 
indicate that a great algebraic simplification is attained if C2 = 0. In fact, in this case we can re-scale the function 
F, as F ^ ^aV(a+b)^ j^f^gj. ^}jich the solution (QT]) with C2 = reduces to 

dS"^ = FVdt^ - APP^dx^ - BPP^dy'^ - CF^^dz^ - ^^^c^V^ (34) 

where pi,p2,P3 are the parameters introduced in (j23p . In sections 3, 4 and 5 we use this line element, which we call 
Kasner-like, for illustrating the fact that physics in AD determines F. 



3 Cosmological models in AD. The STM approach 

The aim of this section is to determine F within the context of induced matter theory. To this end we assume an 
equation of state for the effective matter quantities. Our results show that the light-like Kasner metrics (|22|) can 
be used, or interpreted, as 5-dimensional embeddings for a number of cosmological models in AD that are spatially 
anisotropic extensions of the FRW ones. 

For the Kasner-like metric ([M]) . the components of the effective EMT induced on E^, : = i/jq = constant are 
given by (in what follows we simplify the notation by omitting the bar over F in (I34p and the index in ''*^T)„y) 



a^cF 



(a + 6)2yF2' 
a{(7 + /?)(a + 6)F 



2c(a-/3-7)F| 



2{a + b)^VF^ 



a|(7 + a)(a + 6)F 



F _ V 
F V 



2c(/3-a-7)F| 



2(a-|-6)2yF2 



[{a + (3){a + b)F [{ - |] " Ml ^ " - /3)^} 



2(a + 6)2yF2 



(35) 
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Certainly the specific shape of the EMT depends on the embedding. However, there are a number of relationships, 
between the components of the EMT, which are "embedding-independent" . These are 

(7 - /3)Tl + {a- 7)T| + {13- a)Ti = 0, (36) 

and 

(a + 7)ri - (/? + 7)T| = (/3~a)To", 
(a + /3)Ti - (/? + 7)T| = (7-a)T°, 

(a + /3)r| - (a + 7)T| = (7-/3)^0°. (37) 

Let us notice some particular cases: (i) If two of the parameters are equal to each other (axial symmetry), say 
13, then = T^; (ii) If = but a 7^ /3, then = r| = T| = ~T°; (iii) If a = then T| = -T^; (iv) 
In the case of isotropic expansion (a = /3 = 7) then = T2 — (but not necessarily = T| = Tg = —Tq). 



3.1 Perfect Fluid 

Let us consider the case where the effective EMT behaves like a perfect fluid. From ([35]) we find that = Tf = 
requires 



), (38) 

which implies V cx i5^i^4c/(a+6) Substituting this into ^ and using idV /d%lj) = [dV /dt) we find that F must 
satisfy the equatiorjl 

(a + b)FF' + 4.CFF' = 0, (39) 

from which we get 

F = [/(t) + g{^^)-\i-+b)/i^-+-+^) ^ (40) 

where /(t) and g{^) are arbitrary functions of their arguments. The effective energy density p^'^^ = Tq and pressure 
pieff) ^ ^ _r| = -T| are given by 



Jeff) ^ (eff) g^Q Jeff) ^ 

3.2 Ultra-relativistic matter and radiation 

It is well-known that in the case of radiation as well as for ultra-relativistic matter (i.e., particles with finite rest 
mass moving close to the speed of light) the trace of the EMT vanishes identically. From ([55]) we find that T = 
+ Tl + T| + T| = requires 

(« + ^)(f-^)+2c(f)=0. (42) 

This equation is the analogue of Following the same procedure as above we find F = [f{t) + g(^)]("+^)/(2'=+a+f')_ 
Therefore, the solution for radiation-like matter resembles that of perfect fluid in the sense that the effective stresses 
p'x^^'^ = —Ti, v'v^^^ = ^2^2 ; V'z^^'^ = are proportional to the energy density, viz., 

where nx,ny and are constants satisfying + Uy + — 1. If we average over the three spatial directions, this 
is equivalent to saying that the equation of state is p^'^ff '' = p'-'^^^^S, where p^^ff '' = —Tl/Z. 



^We note that l|26t remains invariant under the re-scaling F — > /(<^+''). 
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3.3 Barotropic linear equation of state 

For the sake of generality, and in order to keep contact with isotropic FRW cosmologies, let us study the scenario 
where the effective mater is barotropic, that is the ratio is constant. Thus we set 

p(e//) ^ ^p(e//)^ n = constant, (44) 

which for n = 1 and n ~ 1/3 gets back the above-discussed perfect fluid and radiation-like scenarios, respectively. 
Substituting into ((35|) we obtain an equation similar to ([38|) and ((42|l . but with the coefficient (1 -|- 3n) in front of 
the term cF/F. Consequently, V cx . The condition (dV/dilj) = (dV'/dt) then requires 

F = [fit) + .g(V')] »+''+(3"+i)<= . (45) 

Substituting this expression into p4p and making the coordinate transformation [dj / dt)dt — > dt, (dg/d^p)d^p dtp, 
the line element in 51? can be written as 

dS^ ^ - AH'-'^dx' - BH^^^dy^ - CH^'''>dz' - (46) 



where H = yi + E-ipj "+''+<^"+^''= E \s an arbitrary constant for n = 1/3, but E = ±1 for any other n ^ 1/3; is a 
positive constant introduced for dimensional considerations. 

3.3.1 Kasner universe in 51? 

We immediately note that the case where E — Q, which requires n — 1/3, gives back the well-known Kasner universe 
in bD. In fact, setting t cx T'^<^+b+-2.c)/(a+b+c) ^j^g ^j^g element dH]) reduces to 

dS^ ^dT^ - At^i^ dx^ - Bt^i^- dy^ - Ct'^'^^ dz^ ± Dt^"^ dip^ , (47) 

where 

aa Pa ja c 

91 = 7 — TT"^' 92 = 7 — , , , 93 = 7 — , , , 94 = , , , : (48) 

(a + + c) (a + b + c) [a + b + c) a + b + c 

satisfy S^^j^g^ = Y^f^^qf — 1, typical of the Kasner universe in 5D. In order to avoid misunderstandings, it may be 
useful to reiterate our terminology: (i) the light- like Kasner metric is (P^ . which depends on the light-like variable 
u; (ii) by Kasner- like metric we refer to (j34[) . which depends on one arbitrary function of t and ip, and (iii) Kasner 
metric is the usual name given to (|47p . which depends only on t. 



In general, for the 4-dimensional interpretation of ()46|) we should notice that on every hypcrsurfacc i/j ^ tpo = 
constant {ip = ipQ = constant) the proper time t is given by 

VDdi (3n + l)c 

dr ^ ± — , TO = TTF ^1 7-^ TT- (49) 

(t-fSi^o)™ 2[a + b+{3n+l)c] ^ ^ 

Bellow we consider several cases. 

3.3.2 Anisotropic Milne universe 

If m = 0, theiH 

n=-i. (50) 



'We exclude c = because it corresponds to empty space, i.e., T^^ = 0. 
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In terms of the proper time r = \/D{t + Eijjo), the metric induced on 4-dimensional hypersurfaces can be written 



as 



ds^^dSl = dr^ - ArP'dx^ - BrP^dy^ -CrP^dz^, (51) 



where pi are the parameters defined in (|23p and A, B and C are some new constants. In addition, rii, the ratios of 
the anisotropic stresses to the energy density (1431) are given by 



(a - /? - 7) (/3 - a - 7) (7 - a - /?) 

(a + fe) 



n, = ^ = ^ SttGp^-^/^) ^ ^"^"^ ^ (52) 



For a = /3 = 7 we find Ux = = = —1/3 and consequently we recover Milne's universe, as expected. 

3.3.3 Anisotropic de Sitter universe 

If m = 1, then 

1 2(a^ + /3^ + 7^) 

3 3(a^ + 07 + /37) ■ ^ ^ 

From we get (t + -B'i/'o) oc e*^/^. Taking the negative sign, the induced metric in AD can be expressed as 

= dSl^^ = dr^ - AeP'^/^dx^ - BeP^^/^dy^ - CeP'^^^^dz^ (54) 

For this metric we find 

nx^-i^^±J^±M, n, = -^J^^±J^±^, n.^J^^±J^^±^, S.Gp^^ff^ = ^ \ . (55) 
c c c (a + o)^Z3 

In the case of isotropic expansion {a ~ [3 ~ ^) these equations yield Ux ^ ny ^ Uz = n = —1 and (|54p 
reduces to the familiar de Sitter metric with cosmological constant A(4) =3/1?. An interesting conclusion here is 
that an anisotropic universe can enter a phase of exponential expansion (inflation), without satisfying the classical 
"false- vacuum" equation p — —p (see (IS^l) ). 



(56) 



3.3.4 Anisotropic power-law FRW universe 

For m 7^ 1, from we obtain 

where tq is a constant of integration. Thus, the induced metric in AD becomes 

ds^ = dSl = dr^ - Ar^^dx^ - Br'^'^dy^ - Cr^'^dz^, (57) 

It/- 

where we have set tq = 0; A, 5, C are constants with the appropriate units, and n is given by 

^ ^ 4« ^ 4(a + /3 + 7) , . 

2(a + 6) + (3n+ l)c 2{a^ + (3^ + -f^) + {3n + l){al3 + aj + (3j)' ^ ' 

We note that the denominator of k is non-zero because by assumption here m ^ 1, see (l53l) . The effective matter 
quantities are 

p(-ff) =np^-ff\ 8nGp^'^f^ = ^ (59) 

4r 

We also find 

2a-)- (3n- l)(/3 + 7) 2/3 (3n - l)(a -t- 7) 27 + (3n - l)(a + /3) 



2a " 2a 2a 



(60) 



We note that for c = 0, the space is empty (p^^-^^-^^ = 0), and the line element ([57)1 yields the well-known Kasner 
solution in 4D. Besides, for n = 1 and n = 1/3 the above expressions reduce to those obtained for perfect fluid and 
radiation-like matter discussed in sections 3.1 and 3.2, respectively. 
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3.3.5 Isotropic expansion: spatially flat FRW universe 



The above expressions evidence the fact that for anisotropic expansion the effective EMT behaves like a perfect fluid 
only for n — I. In contrast, isotropic expansion allows perfect fluid for any value of n. In this case the 5D metric 



(|46)) can be written as (we omit the tilde over t and t/j) 

dS^ - ^^^^ C(t + ^^2/(2+3n) r , 2 , , 2 , , 2] ^^1^'^ /g, N 

For n 7^ — 1, on every hypersurface S^, it reduces to 

ds^ - dT^ - Ct4/3("+i) [dx^ + dy^ + dz^] , (62) 
which is the familiar flat FRW model with perfect fluid 

4 

p = np, 87rGp=— — ^ ^ . (63) 

For n = —1 we recover the de Sitter spacetime as shown in (j54p . 

To finish this section we would like to emphasize that although the metrics with a = and c — correspond to 
empty space (Ricci-flat in 4Z3), they are different in nature. For a — the spacetime is Minkowski (Riemann-flat) in 
5D and 4_D, while for c = the components of the Riemann tensor are nonzero in 5D and in the 4D subspace E^. 



4 Cosmological models in 4D. The braneworld approach 

The preceding discussion shows that, in the framework of STM the Kasner-like metric (p4| embeds a large family of 
AD cosmological models that are anisotropic versions of the FRW ones. However, one could argue that the effective 
matter quantities (|35p do not have to satisfy the regular energy conditions [36j , or any physically motivated equation 
of state, because they involve terms of geometric origir0. 

In this section we will see that the 5D metric (I34p can be completely determined if one imposes an equation of 
state on the matter in the brane. Although the concept is the same as in section 3, the physics here is different. 
Namely, in this approach the spacetime is a singular hypersurface and, for the bD Kasner metrics under consideration, 
there is an effective non- vanishing cosmological term in AD (the brane). As a consequence, the time evolution as 
well as the interpretation of the solutions in AD is distinct from the one obtained, under similar conditions, in the 
framework of STM. 



4.1 The braneworld paradigm 

In order to make the paper self-consistent, and set the notation, we give a brief sketch of the technical details that 
we need in our discussion. In the simplest RS2 braneworld scenario our universe is identified with a fixed singular 
hypersurface "S^^ (called brane) embedded in a 5-dimensional bulk with Z2 symmetry with respect to the brane. 
The discontinuity of the extrinsic curvature across is related to the presence of matter on the brane, which is 
described by an EMT that we denote as T^,y. Thus, now the Einstein field equations in 51? are Gab ~ ^^5)^ab""'^''' 
where fc^g^ is a constant with the appropriate units and t^™'"^^ — S'^6gT^^6{ip)/^. 
Israel's boundary conditions I37j relate the jump of to r^i,, namely, 

kl) 1 

{Kf,^,^+ ~ K^^.^- ) = -e-^{T^^ - -rg^^). (64) 
Now, the assumed Z2 symmetry implies = —K^j_y^^- . Consequently, 

' *6 ^ ^b 

'^In fact, the effective EMT defined by JSJ contains a contribution, given by -Bq^, which is the spacetime projection of the 5D Weyl 
tensor and connects the physics in 4D with the geometry in 5D. 
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"(5) 



(65) 



where the extrinsic curvature K^^i, has to be evaluated on S^^. From G^4 = it follows that r/f.^ — 0. Thus t^^ 
represents the total, vacuum plus matter, conserved energy-momentum tensor on the brane. It is usually separated 
in two parts [38] , 



^9 '^i^ ) 



(66) 



where a is the tension of the brane, which is interpreted as the vacuum energy density, and 7^^ represents the 
energy- momentum tensor of ordinary matter in AD. 
From ([65]) and (|66)) we get 



K., = - 



2- [T^.^^9^.{r + a) 



Substituting this expression into ^ we obtain [12j 
where 



A 



(4) 



, (5)^ 

^ 12 ' 



'- (5)^ 

6 ' 



and 



-T T" 



1 

24' 



(67) 

(68) 

(69) 
(70) 

(71) 



All these four-dimensional quantities have to be evaluated on S^^. They contain two important features; they give a 
working definition of the fundamental quantities A(4) and G and contain higher-dimensional modifications to general 
relativity. Namely, local quadratic energy- momentum corrections via the tensor n^,y, and the nonlocal effects from 
the free gravitational field in the bulk, transmitted by E^v 

4.2 Matter in the brane. Gaussian coordinates 

In the braneworld literature the use of Gaussian coordinates in quite common. In these coordinates the function F 
is given by which under the re-scaling F F°- /('^+'') becomes F — l{t) + Tph{t). 
If we locate the brane at V' = 0, then the metric of the bulk is given by: 

1. For ■(/' > 



-df - A [l{t) + iph{t)f' dx^ - B [l{t) + Tphit)]^'' dy^ ~ C [l{t) + tphit)]"' dz^ - (72) 



2. For V < 



dSl, 



l-iPh 



dt' - A [l{t) - ^h{t)Y^ dx' - B [l{t) - iph{t)r dy' - C [l{t) - ^h{t)Y^ dz' - d^/ . (73) 



Using ^ we calculate the non- vanishing components of K^i, = Kf^^^^+ . These are 



Km = -r^, Kii = — — — , K22 = — — r; — , ^33 = 7 — — Tx — • 

h'' [a + b) (a + 0) (a + 0) 



(74) 
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We assume that the matter in the brane satisfies the equation of state 

P = '>^P^ (75) 

where p = , p — {p^ + Py +Pz)/-^ and p^ ~ ~Ti , Py ~ ^^2^ Pz = ^'^i ■ Using these expressions, from (1551) . with 
e = — 1, and we obtain 



(1 + n) 
2 

(1 + ^) 



h (2 + 3n)a2 fh 
1^ 3(a + 5) \1 



3(a + 6) 



We notice that in cosmological appHcations the metric function g^Q is subjected to the condition 



500 1, 



= 1. 



Thus 



h{t) = s i{t), s = ±1. 



(76) 



(77) 



(78) 



Therefore, we have two equations for the three unknown a,p and l{t). Taking the covariant divergence of (j66p . it 
foUows that to conserve both the total brane energy-momentum tensor r^i^ and the matter energy-momentum tensor 
7^1/, we must have u = ctq = constant. Then, using ([78)1 we integrate the first equation in (l76l) and obtain the scale 
factor a33 



lit) 



V ■ 



3(a + 6) 



(79) 



where Ci and C2 are constants of integration. We note that rj is positive for arbitrary values of a, /3, 7 and n > —1. 
Therefore, if we choose s = — 1 and set C2 = 0, then the "origin" Z = is located at t = —00. Thus, from ((75|) we find 
p cx (T = (To for all t, regardless of the value of n (but n ^ —1). The resulting metric on the brane is de Sitter-like, 
with different rates of exponential expansion in every direction, similar to the models discussed in (|54p . 



4.3 Non-Gaussian embeddings 

The question may arise of whether the simplicity of the above scenario is not a consequence of the simplifying 
assumption of Gaussian coordinates. In order to investigate this question, we consider here the embedding that 
arises from the choicaj 



V(x 



ppqc/{a+b) 
P'2 ■ 



where q is some constant. With this choice the metric in 5D can be written as 

pqc/ (a+b) p2 



pi2 



Now, using [dV/dip) = (dV'/dt) we find 



F =m+i;hit)]'^''+'^/'^''+'-'''\ for 



a + b 



(80) 



(81) 



(82) 



^From I I70I I. with e = —1, it follows that a must be positive in order to ensure G > 0. 

^This is a simple combination between V oc _Fi<'['^('^+3")/(''+'')l for the anisotropic FRW models considered in Section 3.3, and 
V oc F/F''^ for the Gaussian embedding discussed in Section 4.2. 
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and 



F = l{t)e 



i>h{t) 



for q = 



(83) 



where l(t) and h{t) are arbitrary functions of integration. Again, if we locate the brane a,t ip — 0, then the metric in 
the Z2-synimetric bulk is obtained by replacing -0 ^ 1-01 in ([5^ and Following the steps used in Section (4.2) 

we find 



2s 



1 



I (2 + 3n)a2 + 3qc I 
1^ 3{a + b-qc) [l 



(84) 



It is interesting to note that in the case where q — {a + b) /c, the equation for a can formally be obtained from ((84l) 
by setting q = 0. Consequently, ([55)1 yields models on the brane that are identical to those in Gaussian coordinates 
([7S)l . although the metric in the 5D bulk is completely different in both cases. 

The conclusion emanating from (j84p is that, within the context of the 5-dimensional Kasner spacetimes under 
consideration, Gaussian and non-Gaussian embeddings generate the same physics on the brane. In particular, the 
assumption of constant a, which is equivalent to a constant cosmological term A(4), obliges the universe to expand in 
a de Sitter anisotropic form regardless of (the choice of) the embedding. This is quite analogous to the cosmological 
"no-hair" theorem/conjecture of general relativity. 



4.4 Vacuum solutions on the brane 

Since the extra dimension is spacelike, the solutions to the field equations are invariant under the transformation 
(a;, y, z) ?/;. However, the physics in AD crucially depends on how we choose our ordinary 3_D space. 

In order to illustrate this, let us permute ip ^ z in the solution given (|8T|) and ((82| . Also, to avoid misunder- 
standing we change F{t,4') H{t, z). Using this notation, we find that the metric 

dS'^ = 2 - AHP^dx'^ - BHP^dy^ ~ m^'^'^+^^dz'^ ± CHP'dip^, (85) 

^ z 

where Hz = dH/dz and H = [l(t) + zhlt)]''"'^''''^'""'^'' is also a solution of the field equations Rab = 0. Although 
([81]) and ([85]) are diffeomorphic in 5D, their interpretation in 4D is quite different. Specifically, unlike ([81]) in ([85^: 
(i) the extra dimension can be either spacelike or timelike, (ii) the spacetime slices are non-flat, and (iii) the 
metric of the spacetime is independent of ip. As a consequence of the latter, the extrinsic curvature defined by 
([4|), vanishes identically. Which in turn, by virtue of (|65p . implies t^^ — 0, i.e., the spacetime (the brane) is devoid 
of matter (7^^ — 0) and A(4) = 0. 

Clearly, other 51? metrics with properties similar to (|85p can be constructed from the solutions ([^5)) - ([55|) of 
section 2 as well as from ([10)) and ([^5| of section 3. The conclusion here is that the spacetime part of the 5D Kasner- 
like metric (I34p . after the transformation ip <-!■ {x,y,z), can be interpreted as vacuum solutions in a braneworld 
Z2— symmetric scenario. 



5 Static embeddings 

As we noted above, when a 5D metric is independent of ip, the extra dimension can be either spacelike or timelike. 
This is a general feature of the 5D field equations [40]. Therefore, after the transformation tp t the 5D metric 
still satisfies the field equations Rab = 0. The interesting feature here is that after such a transformation the line 
element induced on AD hypersurfaces is static, instead of dynamic as in sections 3 and 4. 

A simple 5D line element that illustrates this feature, in a quite general way, can be obtained from the Kasner-like 
metric p4p in the synchronous reference system. In fact, making the transformations z^ F{t,ip)^H{t, z)\ and 
it ^t, H{t, z) VF(0, z), from ([31]) and ([33|) we obtain 

dS'^ = CWP'^df - AWP'dx^ - BWP^dy"^ - dz^ + c^V^ (86) 
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where (We recall the re-scalling of F introduced at the end of section 2.) 

W = M{z) + tPN{z). (87) 

This is a solution of the 5D equations Rab — for any arbitrary functions AI{z) and N{z). It explicitly depends 
on the extra dimension ip, which now is timelike. At this point it is worthwhile to emphasize that in modern 
noncompactified 5D theories both, spacelike and timelike extra dimensions are physically admissible [9]. 

Once again the choice of the functions M{z) and N{z) depends on the version of 5D relativity we use to evaluate 
the properties of the matter content in AD. Bellow we illustrate this by considering the induced matter approach, 
used in section 3, and the braneworld paradigm used in section 4. 



5.1 Static solutions with planar symmetry in conventional AD general relativity 

It is not difficult to show that the components of the effective EMT, induced on spacetime hypersurfaces : ip = 
■00 — constant, for the metric ([55)1 satisfy algebraic relations similar, but not identical^, to those in ([55]) and (|37p . 
which are independent on the specific choice of M and N. We omit them here and present the case where the 
effective matter quantities satisfy the barotropic linear equation of state (|44p . In such a case we find 



Miz) = CNizr-^,Niz), k^- ^ f^^'V2'''\T\^^^\''i, , (88) 

where C is a constant of integration. Thus, in the 5D solution (jSG]) : W — CN^ + {tp — tpo)N, which implies that 
the metric induced in 4D is independent of the choice of the hypersurface . The effective energy density in AD is 
given by 

C2[(3n+l)(a + /3)+27](a + 6)2' ^"^^ 



and the stresses are 

P'x^^'^ _ (3n - 1)7 - a{in + 1) p^^^^ _ [in - 1)7 - /3(3n + 1) p'f^^'^ _ (3n + l)(a + /3) + 27 
~ 27 ' ~ 27 ' ^(^771 ~ 27 ■ 



(90) 



We note that for fc = 1 these quantities are constants and p^*^-^^' < for all values of a, (i and 7. In what follows we 
assume fc 7^ 1. 

Since 3331^, = ~C^k'^N^''''^^^dN/dzy we can make the coordinate transformation N'-'^^^^dN dz, i.e., N ~ 
In terms of this new coordinate the static metric in 5D is generated by (henceforth we omit the bar over z) 

Wit, ij)^C + {^- V'o)^'/''^-'^ 1 (91) 

The matter quantities induced in 4D decrease as 1/z^. Therefore, the above equations represent static "pancake-like" 
distributions where the matter is concentrated near the plane z — 0, while far from it p —^ 0. 

Except for the singularity at z = 0, the matter distribution presents "reasonable" physical properties. Indeed, 
for every value of n, the "physical" conditions p^'^^^^ > and p > are satisfied in a wide range of parameters 

a, /? and 7. As an illustration, in the case of axial symmetry with respect to z, for n = these conditions hold in the 
range -2/3/3 < 7 < -/3/2 (a = /3 > 0) or -/3/2 < 7 < -2/3/3 (a = /3 < 0). For n = 1/3, they hold if -2/3 < 7 < -/3 
(a = /3 > 0) or — /3 < 7 < —2/3 (a = /3 < 0). A similar analysis can be extended for other values of n. 

A simpler solution can be obtained from the above expressions in the limiting case where fc = 00, which occurs 
for (3ri + l)(a + /3 — 7) + 47 = 0. In this case (|91l) simplifies to W = (7 [z + (-0 — V'o)] and the matter quantities 
are obtained from (|89p . ([50]) by replacing 7 [(3ri + X)(a + l3)/'i{n — 1)]. In the case of axial symmetry, the line 
element becomes independent of the parameters a, (3, 7 and depends only on n. The effective matter quantities satisfy 
pieff) > Q a,nd p > \px,y,z\ for any n in the range — 1/3 < n < 1/3. 



^For example is now replaced by -(a + ■y)T^ + (/3 + 7)T| + (/3 - a)T^ = 0. 
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5.2 Static solutions on the brane 

We now proceed to use the braneworld technique for evaluating the matter quantities. If we locate the brane at 
ip — 0, then the metric in the Z2-symmetric bulk is generated hy W — M{z) + \-4>\N{z). From ([SS]) . with e = 1, 
and ((66|) we obtain the components of 7^,y. Now the barotropic equation of state ((75|) yields a differential equation 
linking M(z), N{z) and a, which can be easily integrated for constant vacuum energy, u = ctq- Namely, we obtain 

3fc?5^C^o(a + + n)M(z) a[(2 + 3„)(c + ,3) + 3^1 ^ ^ 

^' 2(5 + 37i)[72+7(a + /3)]+4(2 + 3n)[a2+/3(a + /3)] ^> ' ^""^^ 

where is a constant of integration. Using this expression we obtain 

klsp = ■ ^ H — = , (93) 

(a + 6)(2 + 3n)M(z)'= k{a + b){2 + 3n) 



(94) 



where 

~ _ (5 + 3n)[-f^ + 7(a + + 2(2 + 3n)[a^ + f3{a + /3)] 

We note that p — constant for 7 = 0. Therefore, in what follows we assume 7 > 0. Since M(z) is an arbitrary 
function, without loss of generality we can choose it a 

M{z) ~ z^/^ (95) 

which is suggested by the decrease of the effective density discussed in section 5.1. It is not difhcult to see that p 
is positive for a large number of values of a, (3,^. The positivity of the first term is guaranteed by the constant of 
integration E. To illustrate the positivity of the second term, we once again consider the case with axial symmetry 
with respect to z. In this case we find 

lim = 2ao(3/3 + 7)(/3-7) .q.^ 

,^±^P 7(5 + 3n)(2/3 + 7) + 6/32(3n + 2)' ^ ' 

which is positive for any /? > 7 and n > —2/3. 

The main conclusion from this section is that regardless of whether we use the braneworld paradigm or the 
induced matter approach, the basic picture in AD is essentially the same. Namely that the AD part of ((86)) represents 
static pancake-like distributions of matter. 



6 Summary 

The vacuum Einstein field equations for the 5D FRW line element (|6]) allow complete integration in a number of 
cases. In particular for $ = 1, or n = 1, the {t, '0)-component of the field equations provides a relation that leads to 
a set of first integrals [15], [l6j. However, for the simplest anisotropic extension of namely, the diagonal Bianchi 
type-I metric 

3 

dS^ = n^{t,i^)dt^ - ^6,(^,1/;) (d2;')%e$2(^,-0)dV>^ (97) 

i=l 

this procedure does not work. (For a discussion, and a new point of view in the context of braneworld, see |41].) 
Here we have pointed out that making the coordinate transformation du cx [ndt — <f>dV')j dv oc [ndt + ^dip), in 
with e = — 1, the field equations allow complete integration in several physical situations, viz., ^U^ . (|A-7p . (IA-8[) . 
The AD interpretation of the 51? solutions requires the introduction of coordinates adapted to spacetime sections . 

We note that in the case under consideration (7 ^ 0), k never vanishes. In fact, for real parameters a and f3 the quantity [a^+/3(cj+/3)] 
is always positive. On the other hand, [7-^ + 7(0 + /3)] = requires o + /3 + 7 = 0, i.e., a = 0, which corresponds to Minkowski space in 
5D and AD. 



16 



We introduced the (t, tp) coordinates by setting u = F{t, tp) and v — V{t, ijj), and used a foliation of the 5D manifold 
such that is a hypersurface of the foliation that is orthogonal to the extra dimension with tangent = 6f/^. 
From a mathematical point of view the functions F and V can be arbitrary, except for the fact that they have to 
satisfy (|26p . However, from a physical point of view, they are related to two important aspects of the construction 
of the spacetime, namely: (i) the choice of coordinates in 5D, e.g., Gaussian normal coordinates adapted to S^, and 
(ii) the formulation of physical conditions on the matter fields in AD, e.g., some an equation of state. 

Our study shows that there is a great freedom for embedding a AD spacetime in an anisotropic 5D cosmological 
model. Similar results but in a distinct context have been found in |41j . To simplify the algebraic expressions, but 
not the physics, in sections 3, 4 and 5 we have devoted our attention to the study of AD spacetimes embedded in the 
light-like Kasner cosmological metric, which is a simplified version of (jST]) . 

We have seen that the simple one- variable line element (|22p can accommodate a great variety of models in AD. 
Indeed, within the context of STM and braneworld theories, we have shown here that the Kasner-like metrics p4p 
may be used or interpreted as embeddings for a large number of cosmological and static spacetimes in AD. Thus, 
apparently "different" astrophysical and cosmological scenarios in AD might just be distinct versions of the same 
physics in 5D |40j . 

This investigation can be extended, or generalized, in different ways. In particular, we have not fully examined 
the possible AD interpretation of the self-similar homothetic solution (pO)) . Neither, have we investigated the solutions 
(jA-7P and (jA-Sp . An important future development here is the question of how these solutions can be applied in 
the generalizations of Mixmaster or Belinskii-Khalatnikov-Lifschits oscillations, as well as other issues mentioned in 
section 1, which appear in theories with one extra dimension. 



Acknowledgments: I wish to thank one of the anonymous referees for a careful reading of this manuscript as well 
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Appendix: Solving the field equations. Part II 

Here we show two more families of analytic solutions to the field equations ([8|) . With this aim we notice that these 
equations are greatly simplified if we introduce the function 



in terms of which (El) become 



AXuv + 
-t- 

Aauv + 



V = e 



XuVy XyVu 



(A-1) 



V V 

V V 



V V 

Adding these equations and using (|A-ip we obtain an equation for V, namely, 

2VV„, - V„V, - 0, 

whose general solution can be written as 



V 



h{u) + g{v) 



(A-2) 



(A-3) 



(A-4) 



where h and g are arbitrary functions of their arguments. Clearly, the self-similar solution discussed in section 2.1 
corresponds to the particular choice 



(A-5) 
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which satisfies (jA-2p and (|A-3|) identically. 

In what follows, as in section 2.1, we introduce a new set of null coordinates u and v by the relations h{u) = ciu 
and g{v) = C2V, where ci and £2 are constants. In terms of these new coordinates V = (ciw + C2v) . Substituting 
this expression into the first of the equations (jA-2p . and dropping the tilde characters, we obtain 

2 (ciU + C2V) \uv + CiAi, + C2\u = 0. 

A similar expression holds for /i. The solutions bellow are obtained under the assumption that e"^ and e'^ are 
separable functions of their arguments. In which case, the above equation implies that they are proportional to 
g±(ci«-c2«)^ The metric function e'^ — V 6^''''+^' automatically satisfies the third equation in (|A-2p and is non- 
separable. Consequently, there are two different families of solutions corresponding to whether cx e^'' or cx e''. 

• In the case where e'*' cx e^^, the field equations Ruu = and — reduce to 

2Au ~ ci (ciM + C2v) A^O, and 2Av - C2 {cm + C2v) A^O. (A-6) 

These equations completely determine the function A and assure the fulfillment of Ruv = 0. The final form of the 
solution is given by 

dS^ = Coe[('=i"+'==^'')'/'*ldMdw - Cie^'^^^^^^^^da;^ - Cae-^"^"-'^^'')^?/^ - (CiCa)"^ (cm + C2vf dz\ (A-7) 

• Following the same steps as above we find that when e'^ cx the solution is 

dS^ = Coe[3('=i"+'=='')'/'*le-^('=i"-"^")dudw - C e('=i"-'=^'') {dx^ + dy^) - C^^ {au + C2vf e-2('=i""'=^")dz2. (A-8) 

In the above line elements (C, Co, Ci, C2) are constants of integration. We note that the resulting solutions are quite 
complicated even in the case where either ci or C2 are set equal to zero. Although this is a great obstacle for the 
analytical interpretation of these metrics in 41?, it allows us to appreciate the simplicity of the self-similar solutions 
discussed in the main text. 
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